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Wind-US Results for the 
AIAA 2 nd Propulsion Aerodynamics Workshop 

This presentation contains Wind-US results presented at the 2 nd Propulsion Aerodynamics Workshop. The 
workshop was organized by the American Institute of Aeronautics and Astronautics, Air Breathing Propulsion 
Systems Integration Technical Committee with the purpose of assessing the accuracy of computational fluid 
dynamics for air breathing propulsion applications. Attendees included representatives from government, industry, 
academia, and commercial software companies. Participants were encouraged to explore and discuss all aspects of 
the simulation process including the effects of mesh type and refinement, solver numerical schemes, and turbulence 
modeling. 

The first set of challenge cases involved computing the thrust and discharge coefficients for a 25° conical nozzle for 
a range of nozzle pressure ratios between 1.4 and 7.0. Participants were also asked to simulate two cases in which 
the 25° conical nozzle was bifurcated by a solid plate, resulting in vortex shedding (NPR=1.6) and shifted plume 
shock (NPR=4.0). 

A second set of nozzle cases involved computing the discharge and thrust coefficients for a convergent dual stream 
nozzle for a range of subsonic nozzle pressure ratios. The workshop committee also compared the plume mixing of 
these cases across various codes and models. 

The final test case was a serpentine inlet diffuser with an outlet to inlet area ratio of 1.52 and an offset of 1.34 times 
the inlet diameter. Boundary layer profiles, wall static pressure, and total pressure at downstream rake locations 
were examined. 
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Workshop Objectives 
(from 2 nd PAW Brochure) 

• Assess the numerical prediction capability (e.g., mesh, numerics, turbulence 
modeling, high-performance computing requirements, and modeling techniques) 
of current-generation CFD technology/codes for Air Breathing Propulsion related 
Aerodynamic flows. 

• Develop practical numerical simulation guidelines for 2-D and 3-D CFD 
prediction of jet related flow fields utilizing Navier-Stokes equations. 

• Explore the underlying physics, flow interaction, jet mixing and dissipation flows 
related to Propulsion Aerodynamics. 

• Enable development of more accurate prediction methods, processes, 
procedures and tools. 

• Enhance CFD prediction capability for practical air breathing propulsion 
aerodynamic design and optimization. 

• Provide an impartial forum for evaluating the effectiveness of existing computer 
codes and modeling techniques. 

• Enhance interests in jet related flows and Identify areas needing additional 
research and development. 

• Cultivate collaboration between the aerospace industry, research institutions 
and academia. 
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Workshop Test Cases 



25° Conical Nozzle 



25° Conical Nozzle 
w/ Splitter Plate 



Dual Stream Nozzle 



_i i i L_ 

15 20 25 30 

X 


Serpentine Inlet (S-Duct) 
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Workshop Format 

• Two-day workshop held at the Cleveland Convention Center 
following the 2014 Propulsion and Energy Conference - July 31- 
August 1, 2014 

• Each group was given 25 minutes to present each of the two 
sets test cases: 

- Nozzle results were presented in Room 6. 

- S-Duct results were presented in Room 7. 

• Organizers will present a summary paper at 2015 Propulsion 
and Energy Conference. 

- Computational results will be consolidated and compared with 
experimental data. 

- Conclusions and observations of trends will be discussed. 

• Select results will be presented at the 2015 AIAA Propulsion 
and Energy Conference. 
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25° Conical Nozzle 


S Instance 1 

- Axisymmetric Conical Nozzle 

• 3.0 in diameter 

• 25° half-angle 

- NPR=1. 4-7.0 (11 cases) 

- Requested data: 

• C</> C v 

• M wa ii on nozzle wall 

• M wa /i from rake in jet plume 

S Instance 2 

- Compare jet plume for NPR=4.0: 

• 25° conical, axisymmetric 

• 25° conical w/Splitter plate 

- Requested data: Flowfield p,T,M,6 

S Instance 3 

- Time-accurate simulation of splitter plate 
vortex shedding for NPR=1.6 

• 25° conical w/Splitter plate 

- Requested data: flowfield snap-shot 
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W nd-US So ver 


• Wind-US v3 

- RANS and Hybrid RANS/LES 
solvers 

- Structured and unstructured grids 

- Numerous turbulence models, 
numerical schemes, and 
boundary conditions 

• All cases: 

- Structured grid solver 

- Roe 2 nd -order physical spatial 
integration scheme (default) 

- Minmod TVD grid flux limiter 
(default) 

- Inflow: p 0 , T 0 held 

- Outflow: p inf held 


• Axisymmetic 

- RANS with SST turbulence model 
(no compressibility corrections) 

- Nominal: CFL#=0.10 on fine grid 

- At= 2.0e-8 s on fine grid for simulations 
showing signs of unsteadiness 

• 3D w/Splitter plate, NPR=1 .6 

- Spalart Detached Eddy Simulation 
(DES) method 

- At= 2.0e-8 s on fine grid 

- DQ limiter on for jet plume zones 

- Zone coupling mode: average 

• 3D w/Splitter plate, NPR=4.0 

- RANS with SST turbulence model 
(no compressibility corrections) 

- CFL#=0.10 on fine grid 

- DQ limiter on for jet plume zones 

- Zone coupling mode: average 
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PAW-Provided Grids 



• Structured 

• Radial grid lines NOT 
perpendicular to centerline in 
plume region 

• (As) wall = 1e-4 inches 

• Zonal boundaries through shear 
layer 

• Zones repartitioned for parallel 
processing 

• Axisymmetric, Conical: 

- 71,424 cells 

- 11 zones 

• 3D w/Splitter Plate: 

- Original grid: symmetry plane aligned 
with splitter plate - NOT favorable for 
vortex-shedding simulation 

- 12,719,400 cells 

- 69 zones 



www.nasa.gov io 



National Aeronautics and Space Administration 


GRC-Developed Grids 



• Created with Pointwise 

• Structured 

• Point-matched zonal interfaces 

• Radial grids lines perpendicular 
to centerline 

• (As) wa n= 1 e-5 inches 

- for y + ~1 at nozzle exit 

• Tighter grid spacing in jet plume: 

- Axisymmetric: Ax<0.05 inches to 
resolve shock and expansion waves 

- 3D: Ax<0.005 inches to resolve 
vortex structures 

• Axisymmetric, w/o Splitter Plate: 

- 352,800 cells 

- 11 zones 

• 3D w/Splitter Plate: 

- 104,112,000 cells 

- 262 zones 
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Convergence 



• Sequence grid: 

- Coarse: every 4 th point 

- Medium: every 2 nd point 

- Fine: all points 

• Axi symmetric, constant 
CFL#: 

- 70,000-120,000 iterations 

• Axisymmetric, constant At 

- 300,000-400,000 iterations 

• 3D, splitter plate, constant 
CFL#: 


Nozzle Exit: C d and C v 



1.02 


0.99 H 


0.98 


Medium Seq. 
Coarse Seq. 


Fine Seq. 


50000 100000 150000 

iterations 


200000 


Jet Plume: Centerline u and TKE 


- 250,000 iterations 

• 3D, splitter plate, constant At 

- DES: 400,000 iterations 


Shown: 25° conical nozzle, 
GRC arid, NPR=2.5. 
Centerline u and TKE took 
longer to converge than C d 
and C v . 



0.03 


0.025 - 



X/D^, 
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Instance 1: Discharge & Thrust Coefficients 




c, 


V 


/-> C r nozzle 1 

2 n \ p-u-r -dr 
P _ Jo 

~ jj > 

P ideal ^ ideal exit 

C ^nozzle I" 2 / \1 

Jd 1 p-« +(p-pj\- r - 

T T nozzle j 

u ue„r J 0 p-u-r -dr 


dr 


Integrated over entire nozzle, including 
base region 

Cd 

- Excellent agreement between two CFD grids 
and experimental data for all NPRs 

Cy 

- Excellent agreement for NPR>3.0 

- Discrepancy between CFD solutions for 
NPR<3.0; no experimental data points for clarity 
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Instance 1: Discharge & Thrust Coefficients 
Comparison with PAW 2012 Results 





• No experimental C v dat for NPR<3.0 

• Compared with PAW 2012 results (71 ,466-cell grid) 

• Still unclear which grid captures C^for NPR<3.0 
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Instance 1: Jet Plume Sonic Lines (Mach=1) 
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0.5 - 


0.4 - 
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♦ Exp: 25° Conical 

♦ Exp: 25° Conical 

♦ Exp: 25° Conical 

♦ Exp: 25° Conical 
Exp: 25° Conical 

Wind-US: Noz25, 

Wind-US: Noz25, 

Wind-US: Noz25, 

Wind-US: Noz25. 

Wind-US: Noz25, 

Wind-US: Noz25, 

Wind-US: Noz25 

Wind-US: Noz25 

-x Wind-US: Noz25, 

-x Wind-US: Noz25, 

-x Wind-US: Noz25 

-x Wind-US: Noz25 

_x Wind-US: Noz25, 

Wind-US: Noz25, 

-x Wind-US: Noz25, 

-x Wind-US: Noz25, 


Nozzle, NPR=2.0 
Nozzle, NPR=2.5 
Nozzle, NPR=3.0 
Nozzle, NPR=4.0 
Nozzle, NPR=5.0 
PAW Grid, NPR=2.0 
PAW Grid, NPR= 2.5 
PAW Grid, NPR= 3.0 
PAW Grid, NPR= 3.5 
PAW Grid, NPR=4.0 
PAW Grid, NPR= 5.0 
PAW Grid, NPR=6.0 
PAW Grid, NPR= 7.0 
GRC Grid, NPR= 2.0 
GRC Grid, NPR= 2.5 
GRC Grid, NPR= 3.0 
GRC Grid, NPR= 3.5 
GRC Grid, NPR=4.0 
GRC Grid, NPR= 5.0 
GRC Grid, NPR= 6.0 
GRC Grid, NPR= 7.0 


' ' i 

0.4 


i ■ ■ 

0.5 


I.'' 

0.6 


x/D : 


jet 


Excellent agreement between solutions of two CFD grids 

Both grids show excellent agreement with experimental data for 
NPR>2.5 
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Instance 2: PAW Grids 



1.5 


0.5 - 


25° Conic Nozzle; No splitter Plate 
PAW Axisymmetric Grid (71,424 cells) 
SST Turbulence Model 
NPR=4.0 
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25° Conic Nozzle; with Splitter Plate 
PAW 3D Grid {12,7 million cells) 
SST Turbulence Model 
NPR=4.0 
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Instance 2: GRC Grids 



1.5 


1 - 


0.5 


25° Conic Nozzle; No splitter Plate 
GRC Axisymmetric Grid (352,800 cells) 
SST Turbulence Model 
NPR=4.0 
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- 1.5 
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-1 - 


25° Conic Nozzle; with Splitter Plate 
GRC 3D Grid (104.1 million cells) 
SST Turbulence Model 
NPR=4.0 


- 1.5 


10 


x/D m 
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Instance 2: Comparison of Wind-US and 
Experimental Shadowgraphs 



NPR=4.0; No Splitter Plate 


NPR=4.0; With Splitter Plate 


25° Conic Nozzle; No splitter Plate 
PAW Axisymmetric Grid (71,424 cells) 
SST Turbulence Model 
NPR=4.0 


25° Conic Nozzle; No splitter Plate 
GRC Axisymmetric Grid (352,800 cells) 
SST Turbulence Model 
NPR=4.0 



25° Conic Nozzle; with Splitter Plate 
PAW 3D Grid (12.7 million cells) 
SST Turbulence Model 
NPR=4.0 


25° Conic Nozzle; with Splitter Plate 
GRC 3D Grid (104.1 million cells) 
SST Turbulence Model 
NPR=4.0 



• Wind-US predicts correct locations of shock and expansion waves. 

• The finer grids do a better job resolving shock and expansion waves. 
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Instance 2: Comparison of Wind-US and 
Experimental Shadowgraphs 



NPR=4.0; No Splitter Plate 


25° Conic Nozzle; No splitter Plate 
PAW Axisymmetric Grid (71,424 cells) 
SST Turbulence Model 
NPR=4.0 



25° Conic Nozzle; No splitter Plate 
GRC Axisymmetric Grid (352,800 cells) 
SST Turbulence Model 
NPR=4.0 






• Wind-US predicts correct locations of shock and expansion waves. 

• The finer grids do a better job resolving shock and expansion waves. 
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Instance 3: Comparison of Wind-US and 
Experimental Shadowgraphs 



25° Conic Nozzle; with Splitter Plate 
GRC 3D Grid (104.1 million cells) 

Spalart DES 
NPR=1.6 



• Wind-US DES indicates a vortex shedding frequency of 28.7-31 .2 kHz, 
somewhat lower than that reported at PAW 2012 (~32 kHz) 

With Medium grid sequence, ~32 kHz vortex shedding frequency 

• No unsteady behavior observed for PAW grid at NPR=1 .6 with DES 

Grid too dissipative? 
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Instance 3: Grid Resolution and Vorticity 



Q 

T3 
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0.04 

0.02 

0 

- 0.02 
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-j 


25° Conic Nozzle; with Splitter Plate 
GRC 3D Grid (104.1 million cells) j 

Spalart DES 
NPR=1.6 


1.05 


1.15 


1.25 


1.35 


X/D 


• GRC Grid solution shown 

- PAW Grid solution: no unsteady flow evident 

• Solutions assume flow is symmetric; only 180° sector modeled. 

• Instantaneous solutions shown. 

• Downstream mesh fine enough to resolving vortex structures. 
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Observations 



• 1 0 of 22 axisymmetric solutions showed unsteadiness when running with 
constant CFL# and required constant time step to converge. 

- Included lower NPRs (<2.5) and high NPRs (5.0, 7.0) 

- Included more GRC Grid solutions than PAW Grid solutions. 

- Convergence required 3-4 times as many iterations. 

• Both grids are good for predicting nozzle performance ( C d and C v ) 

- Some difference between solutions for NPR<3.0, no experimental data for clarity 


• Both grids captured compression and expansion waves 

- Finer GRC grid was superior at resolving and pro-longing waves 

• Curious why the provided 3D grid with splitter plate was configured with the 
symmetry plane aligned with the splitter plate. 

- This seemed a poor choice for observing unsteady vortex formation from splitter plate. 

• Vortex shedding was only observed with DES using GRC grid 

- No evidence of unsteadiness using PAW grid 
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Dual Separate Flow Reference (DSFR) Nozzle 

• Ambient Conditions 

— Pamb~ ^ 4.24 psi 

- T amb =5 20.0° R 

• Fan Stream Inflow 

- Inflow boundary (no profile) with no boundary layer 

- Fan Total Pressure: p 0 ,f af /Pamb varies from 1 .4-2.6 

- Fan Total Temperature: T 0 fan = 530.0° R for all cases 

- Turbulence Intensity = 5%; 

ratio of turbulent to molecular viscosity = 1 .0 

• Core Stream Inflow 

- Core Total Pressure: p 0core = p 0fan l 1 .2 

- Core Total Temperature: T 0 core - 530.0° R for all cases 

- Turbulence Intensity = 5%; 

ratio of turbulent to molecular viscosity = 1 .0 
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Analytical Tools 


• Wind-US v3 

- RANS and Hybrid RANS/LES solvers 

- Structured and unstructured grids 

- Numerous turbulence models, numerical schemes, and boundary conditions 

• All cases: 

- Structured grid solver 

- RANS with SST turbulence model 
(no compressibility corrections) 

- Roe 2 nd -order physical spatial integration scheme (default) 

- Minmod TVD grid flux limiter (default) 

- Inflows: p 0 , 7 0 held constant 

- Outflow: p amb held constant 
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Boundary Conditions 
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Grid 




> 86,259 grid points 

> 11 zones 

> (As) wa n=5e-5 inches 

> (As) exit =5e - 2 inches 


> 2-D Axisym metric, structured, point-matched 

> Created with Pointwise 

> Based on PAW-supplied structured grids 


10 


8 


6 

> 

4 


2 


0 
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Convergence 


• 5 iterations/cycle for all steady-state cases 

• Coarse Mesh: 4,000 iterations 

• Medium Grid: 10,000 iterations 

• Fine Grid: time-accurate 

- 20,000 iterations (if initially run time-accurate) 

- 1 60,000 iterations (if initially run steady-state) 

• Grid zones had to be split and reassembled for convergence 
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Grid Topology 



AIAA 



Would not converge with zones 
split vertically behind nozzle lips 


GRC 
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Fine Mesh Divergence in Steady-State 



Steady-State (CFL#=0.25) 


Time-accurate (Af=5e-9) 
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Fan and Core Stream Discharge Coefficients 



Core Stream 


Fan Stream 




Note: NPR is Fan Stream NPR ( Pojar/Pamb ) 
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Total Thrust Coefficient 
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Fan and Core Stream Total Pressure Rakes 



Core Stream 


Fan Stream 
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Note: NPR is Fan Stream NPR ( Pojar/Pamb ) 
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Fan and Core Stream Static Pressure Taps 



Core Stream 


Fan Stream 


- 9 - 

"Vi 

a 


0.8 

0.6 

0.4 

0.2 

0 


Outer Wall 



18 


19 


18 


19 




NPR=1 .4 
NPR=1 .6 
NPR=1 .8 
NPR=2.0 
NPR=2.2 
N PR-2. 4 
N PR-2. 6 


20 


21 

X [in] 


22 


23 


20 


21 

x [in] 


22 


23 


24 



24 



Note: NPR is Fan Stream NPR ( Pojar/Pamb) 
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Mixing 



Contours: Mach Number 


PO, fctr/Pc 


= 1.8 


Grid is not located where mixing 
layers occur 
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Po,far/Pamb ^ ^ 



Coarse Grid 
Solution 


Fine Grid 
Solution 
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Po,far/Pamb ^ 



Coarse Grid 
Solution 


Fine Grid 
Solution 
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Po,far/Pamb ^ ® 



Coarse Grid 
Solution 


Fine Grid 
Solution 
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Po,far/Pamb ^-0 



Coarse Grid 
Solution 


Fine Grid 
Solution 
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Po,far/Pamb ^-2 



Coarse Grid 
Solution 


Fine Grid 
Solution 
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Po,far/Pamb ^-4 



Coarse Grid 
Solution 


Fine Grid 
Solution 
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Po,far/Pamb ^-6 



Coarse Grid 
Solution 


Fine Grid 
Solution 
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S-Duct Problem Description 



Geometry 

- D, = 133.15 mm (at BL) 

- D 2 = 164.00 mm (at AlP) 

- Area Ratio = 1.52 

- Length = 5.23 * D 1 

- Offset =1.34* D 1 

Flow Conditions 


Boundary Layer 



Total Pressure Rake @ AlP 


- Tested in the R4MA facility at ONERA in 2006. 
Run 1112, Data Point 656 

- p 0 = 88,744 Pa 

- T 0 = 286.2 K 

- massflow = 2.427 kg/s 
(for full 360°) 

- AlP Mach = 0.3549 



Wall Static Pressure Taps 
(cp=0,90,180,270) 

Y 

L. 


v 
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Wind-US 



• Wind-US Version 3 

• Used Grid provided by PAW 

- 20 Million points 

- Split into 385 Zones to increase multiprocessor efficiency and 
reduce computational time. 

- Efficiency: 90 iterations per min (10,000 iterations /1 10 minute) 

• Use symmetry: only model half of the geometry 

• Inflow: Specified p 0 , T 0 , Mach 0.01 

• Outflow: Specified mass flow 2.427 kg/s * 0.5 (symmetry) 

• External Outflow: p 0 12.8712, Mach 0.01 

• Turbulence models 

- Menter Shear Stress Transport (SST) 

- Spalart-Allmaras (S-A) 

• Full description of code features: 

- http://www.qrc.nasa.gov/WWW/wind/index.html 
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Solution Convergence 



• Convergence to “steady-state”. 

- Solutions were shown to be periodically unsteady. 

• Constant CFL of 0.5 

• Ran for 600,000 iterations to eliminate transience and 
show unsteady behavior 

• Last 100,000 iterations was looked at to evaluate 
convergence (plotted Boundary Layer every 10,000 
iterations ) 

• Typical runtime on NASA Pleiades is 124 minutes per 
10,000 iterations (wall time) 
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Solution Convergence on the Fine Grid 

SST 



Boundary Layer 



Wall Static Pressure Taps 




<P 
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Boundary Layer Rake Data 

Boundary Layer 
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0 = 180 ° 
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Streamwise Pressure Variation 

Wall Static Pressure Taps 




0>=0° $=90° <D=180° 
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Circumferential Pressure Variation 




/ 


s/D , =2 s/D , =3 s/D ,=4 
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Symmetry Plane - Fine Grid 




0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 

S-duct at M3A=0.4 
PAW 2014 Grid 
NASA-GRC RANS-SST 



AIP 
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AIP Virtual Total Pressure Rake 



PO Recovery: SST=0.9771 

SA=0.9777 

Note the asymmetry 
in the experiment, 
the computations are 
that of a half section. 


EXP 




*2012 CFD results were provided by Dennis Yoder. 


Mach 




Pi/PiO 


- 1.000 

1 0.995 
0.990 
0.985 
0.980 
0.975 
0.970 
0.965 
_ 0.960 
1 0.955 
T 0.950 
| 0.945 
I 0.940 
I 0.935 
1 0.930 
- 0.925 

y 0.920 
0.915 
0.910 


www.nasa.gov 54 


National Aeronautics and Space Administration 


Observations 

• SST is slightly unsteady while SA is still unsteady but 
to a lesser degree. 

• Incoming boundary layer profiles match well with 
experimental data. 

• There is no significant difference between SA and 
SST predictions when compared with the experiment. 

• 2014 simulations are in good agreement with 2012 
simulations except the 2014 p 0 recovery is slightly 
higher. 

• The distortion values of 2014 and 2012 are not the 
same. 
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